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SUMMARY 

^ Analysis is made of the propulsive thrust and specific impulse 

^ attainable by cooling aircraft surfaces with liquefied or solidified 

gases having vaporization temperatures lower than the equilibrium sur- 
face temperature. It is found that for seme coolants a net thrust is 
obtainable without further heat addition at all flight speeds if the 
vaporization temperature and heat of vaporization are sufficiently low. 

Use of coolant vaporization as the sole means of propulsion yields 
low specific impulses at low speeds. As flight speed Increases, the 
attainable specific Impulse increases; if liquid hydrogen iB used as 
coolant, specific impulses comparable to those of conventional liquid 
rocket propellants are theoretically attainable in the hypersonic speed 
range. This propulsion system provides the possibility of flight with 
cooled aircraft surfaces at extremely high Mach numbers. 

As an auxiliary power source, hydrogen vaporization by aerodynamic 
heating can maintain surface temperatures sufficiently low to preserve 
structural integrity at all flight speeds, and the gas can be ejected 
with a specific impulse in the range of current liquid rocket 
propellants. 


INTRODUCTION 

Rocket propulsion is based on the conversion of liquid or solid 
fuel into gases which are es$elled rearward at the highest possible 
speed. Normally, chemical energy of the fuel is utilized to produce 
the desired gas phase, and to raise the gas temperature to produce 
high Jet velocities . 

In the present report it is proposed that aerodynamic (or atmos- 
pheric) heating be used to vaporize and heat a liquefied or solidified 
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gas having a vaporization temperature below the equilibrium temperature 
of the aircraft surfaces. The resulting gas is to be collected in a 
chamber and expelled through a nozzle as in conventional rocket propul- 
sion. The thrust and specific impulse attainable with such a propul- 
sion system are determined as function of Mach number for a variety of 
possible coolants. Use of coolant vaporization as an independent or aB 
an auxiliary power plant is discussed, and the design parameters that 
determine the efficiency of the system are analyzed. Many engineering 
aspects of the system, however, are not considered herein. 


ANALYSIS 

To arrive at an estimate of the thrust obtainable by the coolant 
vaporization technique, oonsider a missile such as that shown in the 
following sketch. A liquefied or solidified gas having a boiling point or 
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sublimation point considerably below equilibrium surface temperature 
is circulated or stored between the outer and the inner shell. The 
vaporized gases are collected in a chamber C and are allowed to expand 
through the nozzle. The liquid or solid phase is assumed to be at the 
equilibrium vaporization temperature corresponding to the chamber pres- 
sure p c , and the gas phase is assumed to be raised to a .final mean 

temperature equal to the chamber temperature T c . 

The total heat entering the coolant through the external surfaces 
is q^A^ Btu per second. (Symbols are defined in appendix A.) This 

heat vaporizes w pounds per second of the liquefied or solidified gas 
and raises the gas temperature from the vaporization temperature T t 
to a final temperature T 0 . The heat balance is therefore: 


2888 
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^■W^W ” V J^V + °p,C (*0 ” 

where h^ Is heat of vaporization. The mean rate of heat flow into the 
coolant, per square foot of external surface area, can he written 


^w _ kw^e “ *w> - Po ?Q °p,0 St 0 ( T e " ^w) 

where T e = T 0 ^1 + 2Ll_i 6 Mo' 

In this equation, St is the Stanton number, T e is the equilibrium 

surface temperature for an insulated surface, s is the recovery factor, 
and T^ is the mean temperature of the entire cooled surface. 

In terms of the specific impulse I of the gases at chamber con- 
ditions and the total weight flow w, the thrust is 


Thrust = Iw 


(3) 


The net axial force on the missile, in horizontal flight, Is 


F ne t = Iw - Drag = Iw 


(°B,1 + ^ °f) 


p o y q 2a o 

2g 


(4) 


where is the total drag coefficient of the missile (based on Aq) 

with the exception of the friction drag of the cooled surfaces, Cf is 
the mean skin- friction coefficient for the cooled surfaces, and Aq Is 
the maximum frontal area. Use of equations (l) and (2) results in the 
following form of equation (4): 


^net 

I Y 


- °F,net = c f 




(5) 


where G is the ratio of thrust generated per unit area of the cooled 
surface to the friction drag of that unit area, and is expressed as: 


G 


[e£\{ ^ o\ °p,o<yy 

' I £2 Vc f " w VI % j V * «»,.«*. - V 
“ 8 


( 6 ) 


*4 


For turbulent boundary layers, the ratio Si^/i Cj> is found in refer- 
ence 1 to be about 1.2; for laminar boundary layers (ref. 2) this ratio. 
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for Pran&tl number of 0.72, is about 1.245. Inasmuch as these values 
are almost the same, the value 1.2 will be used for both turbulent and 
laminar flow. Similarly, the recovery factor, which is about 0.845 and 
0.88 for laminar and turbulent flow, respectively, will be assumed to 
be 0.88 for both types of boundary layer. With these assumptions, G 
is identical for laminar and turbulent flow. 

Equations (5) and (6) indicate that, for T w < T e , it should be 

possible to reduce the net drag of a missile by utilizing aerodynamio 
heating to vaporize a coolant. In particular, if G can be made 
equal to unity, the friction drag of the cooled surface Is exactly 
cancelled. If G can be made sufficiently greater than unity, it 
should be possible to make the bracket in equation (5) equal to or 
greater than zero, in which case a net thrust can theoretically be 
attained. It Is shown in appendix B that values of G greater than 
unity are compatible with the first law of thermodynamics. 

Evaluation of specific impulse I and thrust parameter G . - The 
quantity gl, which Is the effective jet velocity, can be written 


gl 



(7) 


where m c is the molecular weight of the propellant gas and <p Is 

a dimensionless function of exit pressure p e , chamber pressure p 0 , 

ambient pressure p 0 , and the ratio of speoific heats of the 

propellant y_ : 
c 



where cpmax 



( 8 ) 
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For expansion to ambient pressure (p e = p 0 ), equation (8) reduces 
to 


re - 1 



The maximum value of <P (and consequently the ideal value of I) is 
obtained when p e /p c = 0. Yalues close to this ideal may be attainable 
at extremely high altitudes. This ideal specific impulse is plotted 
in figure 1 as a function of final gas temperature T c , for several 

possible coolants. 

For hydrogen and helium, the specific impulse obtained when 
p e /p 0 = l/6 is also plotted to indicate values corresponding to pres- 
sure ratios that should be readily attainable at lower altitude. 

Of the coolants considered, only hydrogen and helium yield specific 
impulses in the range of current liquid propellants. (Although the 
liquid hydrogen - liquid oxygen propellant systems and a few others 
can attain specific impulses in the range 300 to 400 lb-Bec/lb, most 
current systems are in the range 200 to 300 lb-sec/lb.) The conclu- 
sion to be drawn from figure 1 is that the coolant-vaporization method 
of propulsion can have propellant consumption comparable to that of 
conventional liquid propellants if hydrogen or helium is used and if 
chamber temperatures considerably higher than the vaporization temper- 
atures are attained. 


From equation (6), the following dimensionless thrust parameter 
can be derived: 



= 1.2 



%,o T o \ 

°p, c (Tq ~ Ty)/ 


(9) 


where Tq and ag are ambient static temperature and sound speed, 

respectively. For a given ambient temperature, this parameter depends 
only on the propellant characteristics and an the ratios P e /p c and 

Po/Pe - ■ Etl figu*'© this parameter is plotted against chamber temper- 
ature for the propellants used in figure 1. The ambient temperature 
of the isothermal region of the atmosphere (35,000 to 100,000 ft) has been 
assumed. The data used to plot these curves (as well as those of fig. l) 
are given in the following table, and were obtained from reference 3. These 
values of hy and T 0 correspond to pressures close to atmospheric; 

7 C and Cp j0 are mean values appropriate to the low temperature range. 
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Substance 

h T , 

Btu/lb 

°P,c> 

Btu/lb -°R 

°E 

m c 


A 

67.7 

0.133 

151 

40 

1.67 

CO 

90.8 

.259 

149 

28 

1.41 

C0 2 

160 

.184 

350 

44 

1.3 

He 

10.8 

1.25 

7.2 

4 

1.67 

*2 

194.5 

2.64 

36 

2 

1.5 

^2 

85.7 

.256 

139 

28 

1.4 

°2 

91.7 

.228 

162 

32 

1.4 

nh 3 

588 

.523 

432 

17 

1.31 

h 2 o 

970 

.482 

672 

18 

1.31 


Figure 2 shows that helium, "because of its low heat of vaporiza- 
tion, yields "by far the largest thrust parameter at temperatures near 
the vaporization temperature. At higher temperatures, however, where 
the specific impulse is higher, helium and hydrogen yield considerably 
lower values of the thrust parameter than other substances. These lew- 
values are due primarily to the large specific heats. 

From figure 2, the ratio of thrust to friction drag developed by a 
cooled surface (G) can be computed as a function of Mach number when the 
mean surface temperature and chamber temperature are specified. In 
figure 3(a), results are shown for a chamber temperature and mean sur- 
face temperature equal to the ambient temperature. For this case, 
argon and nitrogen provide the possibility of obtaining a net thrust 
(Gr> l) for Mach numbers greater than about 3.0, while helium and 
hydrogen produce ideal values of & greater than unity for Mach num- 
bers greater than 4.8 and 6.6, respectively. For the lower pressure 
ratio (pq/pc = 1 /6), the Mach numbers for which G equals unity for 

helium and hydrogen are 6.8 and 9.7, respectively. 

Since the thrust parameter of figure 2 either changes little with 
T 0 or reaches a maximum when T 0 = T v , the maximum values of G are 

approximately attained when the chamber temperature and the mean sur- 
face temperature are as low as possible, that is, equal to T t . These 
values, however, are not attainable experimentally, since expansion 
through the nozzle would produce recondensation of the gas. The max- 
imum values of G are nevertheless shown in figure 3(b) as a limiting 
case, and to show the variation of G with Mach number when the mean 
surface temperature T v (assumed equal to T c ) is less than ambient 

temperature Tq. If T Q is expressed in terms of Mach number in 
equation (9), it is evident that G becomes infinite for Mq-»- 0 as 
well as for Mq-*-" if T v < Tq. ■“ 
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For the limiting case shown in figure 3(b), as could he anticipated 
from figure 2, helium produces by far the greatest thrust. Hydrogen, 
argon, and nitrogen, however, all yield maximum values of G greater 
than unity for all Mach numbers if the maximum value of cp can be 
approximately attained. 

Shown in figure 3(c) are values of G obtained when the mean 
surface temperature and the chamber temperature are permitted to riBe 
to 1200° E. Although the specific impulse for this case is much higher 
than for low surface temperatures, the value of G is much smaller at 
each Mach number. 

The conclusion to be drawn from figure 3 is that aerodynamic 
heating can vaporize several liquefied gases at a rate sufficient to 
produce a net thrust at all Mach numbers if sufficiently low values of 
T w and T c are attainable. Such low values of T c , however, produce 

low specific impulse. It is therefore evident that ubo of coolant 
vaporization as an independent propulsion system can be competitive 
with conventional rocket propulsion only at very high Mach numbers, 
where values of G>1 are attainable with higher values of T c . 

Hydrogen or helium vaporization can be used as an auxiliary pro- 
pulsion system (for which large values of G are not required) with 
specific impulses comparable with those attained with conventional 
liquid propellants at all Mach numbers, since the magnitude of T w 
(and consequently T 0 ) is limited only by the requirement that struc- 
tural integrity be maintained. 

Yalues of G and A^/Aq required for propulsion. - Although values 

of G greater than unity were found to be obtainable by vaporization 
of surface coolants at all Mach numbers, equation (5) shows that the 
condition 


^ (G-1)>C DA /C f (10) 

must be satisfied if propulsion by coolant vaporization alone Is to be 
achieved. To estimate the values of G and A^/Aq required for pro- 
pulsion, consider a ccne-cylinder body with cone length-diameter ratio 
fixed at 4.0. The ratio Cp q/C-p for such a body, at zero angle of 

attack, is shown in figure 4 as a function of Mach number for several 
altitudes. (Base drag, boattail drag, and drag due to lift are not 
considered in this example . ) The Reynolds number for these curves is 
based on a length of 50 feet. Friction coefficients for the laminar and 
turbulent cases were obtained from references 2 and 4, respectively. 
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Laminar curves are plotted to one -tenth, the scale of the turbulent 
values. Figure 4 shows that the area ratio required "by condition (10), 
for a given value of G, would be from 7 to 30 times as large (depending 
on altitude) for laminar flow as those required for turbulent flow. 

Shown In figure 5 are the minimum values of G required by condi- 
tion (10) for the case of fully turbulent boundary layer at altitudes of 
35,000 and 100,000 feet. These values are plotted for ratios of wetted 
area to frontal area of 44 and 440. The former value corresponds to a 
cylindrical afterbody of length-diameter ratio of about 9; the latter 
value represents an extreme that might be attained by means of external 
fins. Superimposed on this plot are maximum values of G obtainable 
with hydrogen and helium for several values of the mean surface tem- 
perature T w . In calculating these curves, it was assumed that the 

final gas temperature T c was equal to T v , and that cp = (fi^. The 
maximum specific impulses corresponding to the values of T c are also 
indicated. Figure 5 shows that propulsion at a given Mach number can 
be accomplished with higher specific impulse the larger the area ratio 
and the higher the altitude. For large ratios of surface to frontal 
area, the effect of altitude becomes small. 

These conclusions are further illustrated in figure 6, wherein 
the specific impulses obtainable are cross-plotted against Mach number 
from figure 5. For the large area ratio A^/Aq = 440, hydrogen yields 

higher specific impulse than helium at all Mach numbers. For 
A^/Aq = 44, helium yields higher specific impulses than hydrogen for 

Mach numbers below 4.7 at 100,000 feet and below 6.2 at 35,000 feet. 

In the low supersonic Mach number range, and with A^/Aq = 44, aero- 
dynamic heating alone cannot vaporize liquid hydrogen at a rate 
sufficient to provide a net thrust. 

Figures 5 and 6 illustrate that the specific impulse of the coolant- 
vaporization method of propulsion increases as the ratio of surface to 
frontal area increases, and that the specific impulse increases with 
Mach number. At sufficiently high Mach numbers, the specific impulse 
obtainable with hydrogen vaporization can become comparable with values 
obtained with current liq.uid rochet propellants. 

Heat -transfer requirements . - In the preceding section it was 
assumed that the mean surface temperature T^ was approximately equal 

to the final gas temperature T 0 . Although T w and T c are arbitrary 

parameters in the thrust computation, their aotual values depend on 
details of the internal heat transfer from the surface to the coolant. 



2888 


NACA TN 3140 


9 


From equations (l) and (2), a relation between T w , T 0 , and the 
weight flow w is obtained: 

Ayr P0 c p,0 St 0 (-^e " ^-w) = w + c p,c^c ~ -*vjj (H) 

It is necessary to determine whether a heat exchanger can be devised 
which will produce the desired or assigned values of T w and T c . 

Let A^ be the internal surface area exposed to the coolant (including 
fins and baffles) and let h^ and T^ be the mean heat-transfer 

coefficient from the internal surface to the coolant and the mean 
temperature of the coolant, respectively. Then, since the heat trans- 
ferred into the external surfaces equals the heat transferred from the 
internal surfaces to the coolant, 

A w kwC^e - Tv) = A i hi(T w - Tj[) (12) 

In terms of Stanton number, the ratio of internal to external area 
becomes 

A 1 _ PQ ^0 c p,0 ^ 0(^9 ~ ^y) 

A w p i ^i c p,i ®^i(^w " 

But 

p i T i A Pi = v 

where Ap^ is the internal flow passage area. Substitution of the 

value of w from equation (ll) into equation (14) and the resulting 
value of P^V-l into equation (13) yields 

A i _ + °p, 0(^*0 ~ -^y) 

A Pi " %1 Bt^I, - T ± ) (l5 > 

This equation! shows that the relation between T c and T^ is deter- 
mined by the ratio of internal surface area to internal passage area, 
by the mean coolant temperature, and by the mean internal Stanton 
number. Inasmuch as a portion of the Internal heat transfer takes 
place to a vaporizing liquid or solid, it Is difficult to assign 
appropriate mean values of St^ and T^. However, the order of magni- 
tude of the ratio A^/Ap^ required for specified T w and T 0 can be 

T + T 

estimated by letting T i = -2 —^ — The value of St for heat 


(13) 

(14) 
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transfer to the gas phase in pipes (ref. 5) is roughly 2 to 4xl0“ 3 . If 
St* is assumed to he 3X10"^, equation (15) heoomes 


Api 


333 


\ + °p,0 
G p , c (W - 


(Tq - Ty) 

*0 + T t) 


(16) 


The minimum value of that can he chosen is, as might he expected, 

V. Tq + Ty 

T w 2 - - = Tj. Attainment of this surface temperature requires 

infinite internal surface to passage area ratio. 


The assumption T Q =* T^ yields 


Ai_ 

Api 


667 


c p(^c 


*V 




+ i 


(17) 


For hydrogen, if T c ^ 200° R, the maximum required value of A^/Ap^ 
is about 1000. For T c = 50° R, A±/kp± is about 4000. These values 
are attainable with fins in the annular passage between inner and outer 
shell or by making the height of the annular passage sufficiently small 
relative to its length. It appears, therefore, that the values of T c 
and T w used in figures 5 and 6 should be attainable with suitable 

heat -exchanger design. Pressure drops required in such heat exchangers 
have not been considered. These are functions of the size of the 
missile, the flight Mach number, the altitude, and several other variables. 
At very high Mach numbers, where large coolant mass flows are required, 
the design of a heat exchanger with sufficiently low pressure drop may 
constitute a severe engineering problem. 


Limitations on minimum ratio of nozzle-exit to chamber pressure . - 
Two possible limitations on the minimum attainable value of p e /p c 

exist. One is the requirement that the nozzle-exit pressure cannot fall 
below ambient static pressure without incurring thrust losses. The 
second is the condition that appreciable isentropic expansion below 
equilibrium condensation temperature is probably not attainable. Data 
obtained with air in hypersonic nozzles (ref. 6) indicate that extreme 
purity is required to attain even a few degrees of supersaturation. 

These data also indicate that the pressure-temperature curve for air 
tends to follow the equilibrium condensation curve when the isentropes 
Intersect it. Since the expression given for <P (eq. (8)) assumes 
isentropic expansion, intersection of the isentrope with the condensa- 
tion curve constitutes a limit beyond which I and G cannot be 
calculated by the equations derived herein. 
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To determine the conditions under which the expansion may reach the 
condensation curve "before ambient pressure is attained, the equilibrium 
condensation curre for hydrogen is plotted in figure 7, together with a 
sequence of isentropes. A flight altitude scale corresponding to ambient 
pressure is included on the right of this figure. The second isentrope 
to the right of the condensation curre forms a sort of boundary separating 
ralues of chamber pressure and temperatures for which the altitude limi- 
tation is more significant than the condensation limitation at flight 
altitudes less than about 100,000 feet. For initial chamber temperatures 
and pressures to the right of this isentrope, expansion through the 
nozzle can proceed isentropically to ambient pressure. To the left of 
this isentrope, the condensation curre may be encountered before ambient 
pressure is reached if the altitude is sufficiently high. If, for 
example, the initial chamber pressure is 1 atmosphere (log p = 0) , then 
for altitudes less than 100,000 feet and chamber temperatures greater 
than about 100° B, the ambient pressure is attained before the conden- 
sation curre is intersected. 

To determine the penalty in thrust and specific impulse resulting 
from nonzero ralues of P e /p c and Pq/p q , the ratio of the actual ralue 

of cp to the maximum ralue of <P is plotted for hydrogen in figure 8 
for Po/p e = 1 and 0. These curves were obtained from equations (8) 
and (8a), with y c = 1.5. The ratio p/ y T is also the ratio of actual 

to ideal ralues of G and I, since both these parameters are directly 
proportional to <p. 

It appears from figure 8 that ralues of <P /<P m ax greater than 0.90 
are attainable only with pressure ratios p e /p c less than about 0.007 
if the nozzle-exit pressure is equal to ambient static pressure 
(Po/Pe - l) • ^ "k* 16 maximum chamber pressure is considered to be 

10 atmospheres, ratios of p e /p 0 less than 0.007 are attainable only 

at altitudes greater than 60,000 feet. The numerical conclusions 
reached from figures 5 and 6 must therefore be modified according to 
figure 8 for lower altitudes or lower chamber pressures. 

Effect of specific gravity of propellant on range . - Although 
hydrogen yields the best specific impulse of the coolants considered, 
the density of liquid hydrogen is very low. Hence, a given weight of 
hydrogen will occupy a much larger portion of a given missile volume 
than conventional liquid propellants . To obtain an indication of the 
relative effect of propellant density, consider the range equation for 
flight at constant velocity and constant lift -drag ratio (Breguet path): 



S = T n — I In 
0 D 


( 18 ) 
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where Wo,p is initial propellant weight and We is initial total 
weight less Wq p (empty weight). Let pp he the density of the propel- 
lant, Yp the initial volume occupied by the propellant, and v the 
total missile volume. Then equation (18) becomes 

- I In (1+kPp) (19) 


where 


k = 


v T p T p 


W e 


W P 


For given L/D, V 0 , k, and I, the range is therefore proportional to 
log (l+kp p ). 


If is defined as a volume specific impulse I vo -> , and if 

h Vq 

is expressed in terms of specific gravity, equation (18) becomes 


Ivol - ^ ^ = I In (1+kPp) = I In £l + 62.4 k(specific gravity)j (19a) 


This volume specific impulse, which determines the range at constant 
velocity, constant lift-drag ratio, and constant ratio of propellant 
volume to missile empty weight, is plotted in figure 9(a) for 
k = 2,0 cubic feet per pound and in figure 9(b) for k = 0.10 cubic 
feet per pound. Since range increases with k, large values of k 
correspond to long-range missiles, while small values of k suffice 
for short-range missiles. It Is evident from figure 9 that the range 
for a given propellant volume is very sensitive to specific gravity, 
although specific gravity tends to become less important for long- 
range missiles. For k = 2 cubio feet per pound, liquid hydrogen 
(specific gravity » 0.07) requires a specific impulse of 400 pound- 
seconds per pound to yield a range equal to that of an equal volume of 
propellant with specific gravity 0.7 and specific impulse of 200 pound- 
seconds per pound. These figures indicate that hydrogen vaporization, 
although it provides the best specific impulse, may be less satisfactory 
than other possible coolants If the structural weight and missile size 
required to hold the larger volume are considered. These conclusions, 
based on constant l/d, could be altered considerably if difference in 
fuel weight appreciably affects the optimum L/i) or structural weight 
of the missile . 

Design considerations for optimum range . - If the centrifugal force 
acting on a missile in steady horizontal flight is considered, the range 
equation (eq. 0 - 8 )) becomes (if the rotation of the earth Is negleoted) 
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S = 


TqIL/D 



In 



( 20 ) 


where Kq is the radius of the earth. The lift -drag ratio of a wing- 
body combination can be written 


L _ C Id A d 

35 (°D,d +2c f ,d) A d + °D,b A 0 + ^b^b 


( 21 ) 


where the subscript d refers to wing values and subscript b refers 
to body values. It is assumed that interference between wing and body 
coefficients is negligible and that the body is at zero angle of attack 
regardless of the angle of attack of the wings. Each coefficient is 
based on the area by which it is multiplied, and Cp ^ and Op p refer 

only to drag other than the friction drag. 

For steady flight, the thrust equation (eq.. (5)) becomes 

£ ( 0 - 1 ) - ( 22 ) 

A 0 Of 


where Cp^p includes, as before, all drag except the friction drag of 

the cooled surfaces. If only the body surface Is cooled by the coolant, 
then the friction drag of the wings, as well as their pressure drag, 
must be included in Cp^ If wings and body are cooled, Cp^p is the 

pressure drag of the body and wings. Equation (22) therefore becomes 


Cf 


CD,b 



(23) 


where 

Cp, 2 = Cp^a. + 2C f,d ff the body is cooled 

Cp ,2 = Cp^a if body and wings are cooled 

C f = c f,b A w,b/ A 0 only the body Is cooled 

Cf = Cf £ A^-jj/Aq + 2C fj(i A^/Aq if body and wings are cooled 
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Substitution of equation (23) into equation (21) yields 

L _ ( x ~ g) °L 

D “ Ao 

°D,2 + °D,b 


(24) 


For a given velocity, prescribed aerodynamic coefficients, and weight 
ratio, the range for a coolant -propelled missile is therefore propor- 
tional to l(l-l/G) (eq.s. (20) and (24)). From equations (6) and (7), 
this quantity can be written as follows: 


I - ■= = 


a 0 

g 


<p A fj£3 


M 0 

1.2 


( T T \ 

TO " 




c p,0T0 I 1+0.176 Mq 


(25) 


Plots of (I-l/G) as functions of T c for a given M 0 yield peak 
values of (I-l/G) at an optimum value of T c , to be denoted by 
T c opt . To obtain T c as a function of Mach number, equation (25) 

was differentiated with respect to T 0 . The resulting expression for 

T c,opt ia 


(l+0.176Mo 2 - ^ 
^c ,opt _ \ i 0 / 

T 0 b 2 MQ 2 


(26) 


where 

b 2 = , f 

1.44 <p 2 mQ V C P,0/ 


For a given coolant, the optimum ratio of chamber temperature to 
ambient static temperature is then a function only of Mq and of the 

mean surface temperature T v . Since (I-l/G) increases as T w 
decreases, it is desirable to make T w as small as possible. As 
pointed out in a previous section, a plausible minimum value ifor T w 
is T w = T c . With this value, the egression for T c ^ 0 p^. becomes 

= 1 + (o.l76 + ^ Mq 2 - ~ bM 0 V 4 + (0. 704+b 2 )M 0 2 (27) 


2888 
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For hydrogen (b 2 = 5.364 for <p = nwr) . the variation of T 0 ^ 0pt "with 
Mach number is shewn in figure 10 for an ambient temperature of 392.4° E. 
It is noteworthy that the optimum chamber and surface temperatures are in 
the structurally acceptable range even at Mach numbers of the order 
of 20. 

The maximu m values of (i-l/G) corresponding to these optimum 
temperatures are also plotted in figure 10. With these values, the 
maximum values of I L/D can be obtained from equation (24) when the 
aerodynamic coefficients are specified. With I l/d determined, the 
range can be calculated as a function of velocity and fuel weight from 
equation (20) . 

In evaluating the range obtainable with coolant vaporization 
relative to other systems, it is necessary to consider that the wetted 
area required to yield a net thrust may be larger for the coolant- 
vaporization method than for conventional methods, larger wetted areas 
may have a considerable effect in reducing the maximum lift-drag ratio, 
as can be seen from equation (21) . The effect of this wetted area on 
L/D is included in equation (24). If the resulting l/D is less than 
for a conventional propulsion system, the specific impulse for the 
coolant-vaporization method must be correspondingly higher to yield the 
same range. 


DISCUSSION 

From the preceding sections, it appears that the vaporization of 
liquefied or solidified gases by aerodynamic heating can provide thrust 
greater than the friction drag of the cooled surfaces at all flight 
speeds if the coolant, the heat exchanger, and the pressure ratio are 
appropriately selected. This thrust is attained with no combustion, 
inlets, or machinery (other than, perhaps, a pump) and more than ade- 
quate surface cooling Is provided at all flight speeds. The specific 
impulse associated with this thrust, however, is lower than that 
attainable with conventional rochet fuels unless liquid hydrogen is 
used as coolant and unless the flight Mach number Is very high. 

Detailed investigation of the conditions, or types of mission, for 
which coolant vaporization might be advantageously used as a sole means 
of propulsion Is beyond the scope of the present report. The preliminary 
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analysis presented herein does indicate, however, that flight in the 
atmosphere at extremely high Mach numbers ( >10) can probably be accom- 
plished efficiently with the coolant-vaporization system. Whether Buch 
flight paths are of interest remains to be established. 

Although much of the analysis was concerned with use of coolant 
vaporization as an independent means of propulsion, utilization of this 
thrust source as an auxiliary propulsion system may have more Immediate 
interest. As shown in figure 5, the specific impulse can remain fairly 
high for arbitrarily low Mach numbers if the value of G is reduced as 
the Mach number decreases. Suppose, for example, that the internal 
heat-transfer system is designed to make T w = T c = 392.4° B. Equa- 
tion (17) indicates that such a heat-transfer system (fixed A^/Ap^ 

can maintain this temperature at all Mach numbers if the mass flow w 
is adjusted according to equation (ll) . The specific impulse remains 
fixed at about 257 (fig. l) at all speeds, and the parameter G varies 
as indicated by the straight line T c = 392.4° E in figure 5. The value 

of G remains less than unity for Mach numbers less than about 7. Eor 
a given missile, therefore, the hydrogen-vaporization technique provides 
more and more thrust as the missile accelerates until, at a Mach number 
near 7, all friction drag is effectively cancelled. In addition, the 
surface remains cooled at no cost in propellant weight efficiency if the 
specific impulse of the main power plant is Iobb than 257. Eor higher 
surface temperatures, the specific Impulse obtained with liquid hydrogen 
can be even higher ( see fig . 5 ) . 

An interesting distinction appears from equation (5) between design 
parameters desirable for G>1 and those desirable for G<1. If 
coolant vaporization is to be used primarily as an independent means 
of propulsion (G>l), large values of are desirable so that a 

given net thrust can be attained with smaller values of A^/A q. If 

coolant vaporization is used primarily as an auxiliary power source 
(G<l), then small values of are desirable. For this case, there- 

fore, the possibility that the coolant may be able to stabilize a 
lam i nar boundary layer becomes significant. The curves presented in 
reference 7 indicate that complete stabilization of the laminar boundary 
layer may be feasible in the Mach number range 1.5 to 9 if the surface 
temperature can be reduced almost to the ambient static temperature. 
Although complete experimental verification of this result has not yet 
been obtained, experiments reported in references 8, 9, and 10 indicate 
that the transition Beynolds number is increased as the surface is 
progressively cooled. An evaluation of the relative efficiency of 
cooling with liquid hydrogen rather than with other systems should 
therefore include, in addition to the available thrust, the possibility 
that the friction drag may be reduced from the turbulent to the laminar 
value. 
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For aircraft designed for flight in the atmosphere with a conven- 
tional rocket power plant it appears that hydrogen vaporization can 
provide an extremely efficient cooling system since, as pointed out in 
the previous example, the specific impulse of the coolant flow can "be 
of the same magnitude as that of the main power plant at all flight 
speeds. For use with air-breathing engines, however, it appears 
probable that other more efficient coolant systems can be devised. 

Since the specific impulse for air-breathing engines is much higher 
than for rockets, it would appear to be inefficient to carry large 
quantities of coolant from which specific impulses much lower than 
those of the main power plant are obtained . 

Figure 8 shows that high ratios of chamber pressure to ambient 
pressure are required to attain values of specific impulse close to 
the maximum values. Since high chamber pressures involve heavy 
structures, the high ratios of chamber pressure to ambient pressure 
3 are best attained by flight at very high altitudes. Although the 

computations of the present report were limited, for convenience, 
to the isothermal region of the atmosphere, flight at greater alti- 
tudes by means of coolant vaporization may be feasible and desirable. 
The higher ambient temperatures in the altitude range from 150,000 to 

200.000 feet make this range particularly desirable for coolant - 

. vaporization propulsion. Equation (9) shows that, for a given T c and 

T w , the thrust increases as the square root of the ambient temperature. 
This equation assumes, however, that the skin friction and heat trans- 
fer are related according to continuum flow theory, and that the heat 
transfer from the surfaces due to radiation is negligible. Both of 
these assumptions must be reexamined for altitudes greater than 

100.000 feet. The reduction in skin friction due to Blip flow is 
beneficial if coolant vaporization is used as an auxiliary power 
source, but is hamful if coolant vaporization is to be used as an 
independent power source. Heat loss due to radiation reduces the 
thrust attainable because the mass flow for a given cooled surface 
area is reduced. Eadiation loss can, of course, be kept small even 
at very high altitudes by using surfaces with very low emissivity. 

It is apparent that further analysis is required to determine the 
best flight altitude for utilization of coolant vaporization to provide 
thrust. 


CONCLDBXNG HEMABES 

The preceding analysis indicates that vaporization of aircraft 
surface coolants by aerodynamic heating can, under certain conditions, 
produce a net thrust at all flight speeds . Liquid hydrogen, because of 
its low molecular weight, provides the highest specific impulse at the 
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relatively low chamber temperatures characteristic of this propulsion 
method. Use of hydrogen vaporization as an independent propulsion system 
yields low specific impulses at low speeds; the specific impulse increases, 
however, as Mach number increases and can theoretically equal of exceed 
values obtainable with the best rochet propellants. The specific impulse 
of the coolant-vaporization propulsion system increases with altitude and 
with the ratio of cooled surface area to frontal area. A missile 
designed for efficient use of coolant-vaporization propulsion may there- 
fore be one with thin external fins running longitudinally along the 
body. 


Use of hydrogen vaporization as an auxiliary propulsion system and 
as a means of main taining constant surface temperature shows promise in 
connection with conventional rocket propulsion. Surface temperatures 
sufficiently low to preserve structural integrity can be maintained at 
all flight speeds, and the resulting gas can be ejected with a specific 
impulse comparable to that of conventional rocket propellants. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, May 6, 1953 
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APPENDIX A 


co 
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CD 


Ai 

Aw 

A Pl 

a 0 


c D,b 

^d 

°D,± 

C f 

c E,net 

°L 

C P 

D 

■®net 

G 

g 

h 

h* 


SYMBOLS 

maximum cross-sectional area of aircraft, ft 2 

r 

internal surface area in contact with, coolant, ft* 
external surface area cooled by coolant, ft 2 
internal passage area for coolant, ft 2 
ambient speed of sound, ft/sec 


2 V7q m c/ m 0 ( C V,c\ 

1.2 cp \C p oj 

body pressure drag coefficient 
wing pressure drag coefficient 

total drag coefficient of aircraft lesB friction drag 
coefficient of cooled surfaces 

mean external skin-friction coefficient 
net thrust coefficient 
lift coefficient 

specific heat at constant pressure, Btu/lb-°E 
drag, lb 

net axial force on miBsile 

thrust parameter, defined by equation (6) 
acceleration of gravity, 32.16 ft/sec 2 
heat-transfer coefficient, Btu/sec-ft 2 -°E 
internal mean heat-transfer coefficient, Btu/sec-ft 2 -°E 
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h v heat of vaporization, Btu/lb 

external mean heat-transfer coefficient, Btu/sec-ft 2 -?R 
I specific impulse, lb -sec/lb 

I T0 1 volume specific impulse, sec 

k T p/We 

L lift, lb 

Mq flight Mach number 

m c molecular weight of coolant 

mo molecular weight of air 


P 

P 

E 

*0 

St 

T 

T e 

T t 

t 

Y 


power 

static pressure, lb/ft 2 

mean heat -transfer rate from external boundary layer to 
external surface, Btu/sec-ft 2 

a ft 2 lb 

universal gae constant, 4.97x10 — 

sec 2 lb -mole- °E 

radius of earth, 20 . 9X10® ft 
Stanton, number, = h/pYCp 
static temperature, °E 

T 0 (l + 2L f 1 e Mq 2 ) 

vaporization temperature of ooolant at pressure Pqj 
wing thickness ratio 
velocity, ft/sec 
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P 

W e 

Wo 


W, 


0,p 


w 

X 


a 

P 

7 

e 

P 

<l> 


9max 


total missile volume, ft 3 
initial propellant volume, ft 3 

empty we ight , It 

total initial missile weight, lb 
total initial propellant weight, lb 

coolant flow, Ib/seo 
range, ft or miles 
angle of attack, radians 

Vmo 2 - 1 

ratio of specific heats 
recovery factor 
density, lb/ft 3 

function of 7 C , p e /p c , and Po/p e defined by equation (8) 



Subscripts: 

b pertaining to body 

c values in collection chamber 

d pertaining to wing 

e nozzle exit 

i mean values for coolant 

j jet 

opt optimum 

p propellant 

v vaporization 
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v mean values at external surfaces 

0 ambient stream conditions 
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APPENDIX B 
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RELATION BETWEEN PROPULSIVE POWER AND RATE OP HEAT ADDITION 

The power considerations for the coolant -propelled system are iden- 
tical with those of the conventional rocket engine- 

For steady flight at velocity Vo the propulsive power for either 
system is 


wV.Vq 

Pj = thrust x V 0 = IvVq = — — 


(Bl) 


The portion of the thermal energy of the propellant converted to kinetic 
energy in passing through the nozzle is 

V-2 

=T - JCp,c< T c - V - 4- <B 2 > 

where tj is the propellant temperature following expansion through the 

nozzle. Frcm equations (Bl) and (B2), the ratio of propulsive power to 
thermal power (defined as weight flow times thermal energy used) is 



Equation (B3) shows that, for rocket-type jets, the propulsive power can 
be less than, equal to, or greater than the thermal power of the propel- 
lant, depending on the ratio of missile velocity to effective jet veloc- 
ity. This result Is pertinent to the discussion in the text, since it 
can readily be shown that the rate of energy addition to the coolant 
due to aerodynamic heating Is not sufficient to overcame the friction 
drag of the cooled surfaces at high speeds. The result obtained in 
equation (B3) indicates that, in spite of this inadequacy of the energy 
obtainable frcm aerodynamic heating, values of G greater than unity 
are compatible with the first law of thermodynamics . 

To evaluate the source of the energy corresponding to the propul- 
sive power, the kinetic energy acquired by the propel 1 ant during accel- 
eration to velocity Vq must be considered. The propulsive power does 

work relative to a stationary reference frame; consequently, the total 
energy change of the propellant relative to that reference frame is of 
significance. The total energy per pound of propellant in the chamber 
before ejection Is 
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E = - 2T + *p,d*c < B4 > 

After ejection, the total energy of the propellant is 

(V 0 - V^) 2 . . 

Ej = ' ° ~ 2g J + Jcp^tj (B5) 

The difference between these energies, multiplied by the weight flow, 
is, as might be expected, the propulsive power. Thus, 

wVqV! 

w(E c - Ej) = — — 1 1 * wIV 0 (B6) 

s 

Equations (B5) and (B6) show that part of the propulsive power, 
relative to the stationary reference frame, arises from the kinetic 
energy acquired by the propellant, relative to that reference frame, 
during acceleration to velocity Vq. 

At very low speeds, the kinetic energy of the propellant, aB well 
as the aerodynamic heating and the friction drag, becomes negligible. 

The coolant-vaporization system then develops thrust only if the 
vaporization temperature is less than ambient air temperature. The 
atmosphere then serves as a heat source to vaporize the coolant, 
which develops pressure in the chamber and generates thrust in 
expanding through the nozzle. 
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(a) T w = T c - T 0 . 

Figure 3. - Ratio of generated thrust to friction drag as fraction of 
Mach number for cooled external surfaces. Tg = 392.4° R. 
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Figure 3. - Continued. Eatio of generated thruBt to friction drag as function 
of Mach number for cooled external surfaces . Tq = 392.4° B. 
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Mach number j Mq 
(e) T c = T w = 1200° R. 

Figure 3. - Concluded. Ratio of generated thrust to friction drag as function of Mach 
number for cooled external surfaces. Tq = 392.4° R. 
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Propellant specific gravity 
(to) K = 0.10 cubic .f egt ,per_ pound. 

Figure 9. - Volume specific impulse as function of specific impulse and specific gravity. 
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(a) Optimum chamber (and surface) temperature. 



(b) Max im um value of I - i/G. 

Figure 10. - Maximum value of range parameter (I - l/G) mav and corres- 
ponding values of chamber temperature for hydrogen cooling. T c = T w . 
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